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Abstract. We show that eigenvalues and eigenfunctions of the Laplace-Beltrami 
operator on a Riemannian manifold are approximated by eigenvalues and 
eigenvectors of a (suitably weighted) graph Laplace operator of a proximity 
graph on an epsilon-nct. 



1. Introduction 

Traditionally, discretization in Riemannian geometry is associated with triangu- 
lations and other polyhedral approximations. This approach works perfectly well 
in dimension two but meets a number of obstacles in higher dimensions. It is now 
clear, due to works of Cheeger, Petrunin, Panov and others (see [SJ [T3] H21 [UJ ) that 
in dimensions beyond three polyhedral structures are too rigid to serve as discrete 
models of Riemannian spaces with curvature bounds. In some applications, we get 
a Riemannian manifold as a cloud of points with approximate distances between 
them, see e.g. [TO], [4]. It appears that point clouds arising as discretizations of 
Riemannian manifolds can be effectively distinguished from arbitrary ones. This 
issue will be addressed elsewhere. For triangulations, even the problem of determin- 
ing whether a given simplicial complex is a topological manifold is algorithmically 
undecidable (see e.g. [UJ §6.2] for a simple proof). 

In a few papers, we will try to discuss approximating Riemannian manifolds by 
graphs, of course with additional structures and various boundedness conditions. 
Here we show that the spectrum of a suitable graph Laplacian gives a reasonable 
approximation to the spectrum of the Riemannian Laplace-Beltrami operator. The 
key difference with Unite element and similar methods (see e.g. [7], [5] and an 
interesting recent work [3]) is that in our construction the set of vertices is an 
arbitrary net as long as it is dense enough. There are no local regularity constraints 
and we use only very rough data. 

Let M n be a compact Riemannian manifold (without boundary) and X C M a 
finite e-net. The geodesic distance between x, y € M is denoted by d(x, y) or simply 
\xy\. Given such X and p > 0, one constructs a proximity graph T = T(X, p): the set 
of vertices of the graph is X , and two vertices are connected by an edge if and only 
if d(x, y) < p. In our set-up, we assume that e -C p and p is sufficiently small so that 
p-balls in M are (bi-Lipschitz) close to Euclidean. In addition, we assign weights 
to vertices and edges of T as explained below. Then there is a graph Laplacian 
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operator associated with this structure, see (|1.3[) . Our goal is to approximate the 
eigenvalues Afe(M) the Laplace-Beltrami operator on M by eigenvalues Afe(F) of 
the graph Laplacian. 

This kind of problems were studied before. Fujiwara [9] showed that, if X is an 
e-separated e-net and p = 5e, then the eigenvalues of (unweighted) graph Laplacian 
of the proximity graph after proper normalization satisfy 

C-'XkiM) < A fc (r) < C n \ k {M) 

where C n > is a constant depending only on n = dim M. Belkin and Niyogi [1] 
considered random, uniformly distributed nets in M and showed that, for a suitable 
choice of edge weights (depending on distances), the spectrum of the resulting graph 
Laplacian converges to the spectrum of M in the probability sense. 

In this paper we present a construction that works for an arbitrary net. The 
"density" of the net may vary from one region to another. To compensate for 
this, we need to introduce weights on vertices. These weights determine a discrete 
measure on X and we essentially require that X approximates M as a metric 
measure space. 

The construction. Let e > and X = {xi]f =l be a finite e-net in M . We denote 
by B r (x) the closed metric ball of radius r centered at x £ M . We assume that X 
is equipped with a discrete measure p — ^Pi5 Xi which approximates the volume 
of M in the following sense. 

Definition 1.1. A measure p on X is an e- approximation of volume vol on M if 
there exist a partition of M into measurable subsets Vi, i = 1,...,N, such that 
Vi C B £ (xi) and vol(V,) = pi for every i. 

In this case we also say that the pair {X, p) e- approximates (M, vol). 

Every e-net X in M can be equipped with such a measure. For example, let 
{Vi} be the Voronoi decomposition of M with respect to X and pi — vol(Vi). We 
discuss other constructions and some properties of Definition 11.11 in Section [8] In 
particular, we show that this definition is naturally related to weak convergence of 
measures (see Remark 18.41) . 

Consider the space L 2 (X) = L 2 (X,/j,), that is the A^-dimensional space of func- 
tions from X to K equipped with the following inner product: 

(1.1) (u, v) = (u, v) L 2 {x) = ^ Viu{x % )v{xi), 
or, equivalently, with a Euclidean norm given by 

(1.2) iMpHMii^^^HxOI 2 . 

We think of L 2 (X) as a finite-dimensional approximation to L 2 (M). For the sake 
of brevity, we omit the index L 2 (X) in most formulae in the paper. 

We define the following weighted graph L = T(X, /i, p). The set of vertices is X, 
two vertices x, y e X are connected by an edge if and only if d(x, y) < p. We write 
x y for x, y G X if they are connected by an edge. Both vertices and edges are 
weighted. The weight of a vertex Xi is pi. To an edge e%j = (xi, Xj) we associate a 
weight w(eij) = Wij given by 

2(n + 2) 

W *i = y n n+2 

where v n is the volume of the unit ball in W 1 . Note that wa — w-a. 
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We approximate the Riemannian Laplace-Beltrami operator A = Am by the 
weighted graph Laplacian Ar : L 2 (X) — > L 2 (X) defined by 



The motivation behind this formula is the following. If u is a discretization of a 
smooth function / : M — > K., then the latter sum is the discretization of an integral 
over the ball B p (xi): 



and the normalization constant is chosen so that the normalized integral approaches 
Af(xi) as p — > 0, see Section 12.31 It follows that the graph Laplacian of the 
discretization of / approximates A/ if e <C p <C 1. 

Remark 1.2. One can introduce weights on edges depending on their lengths. For 
example, the above value of uiy could be multiplied by ip(p~ 1 d(xi,Xj)) where ip 
is a nonnegative non-increasing function on [0, 1]. With a suitably adjusted nor- 
malization constant, everything generalizes to this set-up in a straightforward way. 
Probably a smart choice of ip can improve convergence rate. 

The operator Ar is self-adjoint with respect to the inner product ([l.ip on L 2 (X) 
and nonpositive definite, see Section [2Tl Let = Ai(T) < A2(r) < • • • < \n(T) be 
the eigenvalues of — Ar and = Ai(M) < A2(M) < . . . the eigenvalues of — Am- 

Statement of the results. Let Ai = A4 n (K, D, io) be the class of n-dimensional 
Riemannian manifolds with absolute value of sectional curvature bounded by K, 
diameter bounded by D and injectivity radius bounded below by i . 

Throughout the paper, we denote by C various absolute constants whose precise 
value may vary from one occurrence to another (even within one formula) . We write 
C„, Cm, etc, to denote constants depending only on the respective parameters. 

In some of the arguments we denote by er, a\ , . . . , various "small" quantities 
depending on e, p, etc. These notations are local and redefined in each proof where 
they are used. 

The main result of the paper is the following 

Theorem 1. For every integer n > 1 there exist positive constants C' n and c n such 
that the following holds. Let M € M. = Ai n (K, D, io) and T — T(X,p 1 p) be a 
weighted graph defined as above, where (X,(j,) s -approximates {M, vol), p < io/2, 
Kp 2 < c n and e/p < min{l/n, 1/3}. 

Then for every k G Z + such that pAfc(M) < c n one has 



(1.3) 





Hj{u{xj)-u{xi)). 




(1.4) |A fe (T) - \ k (M)\ < C n {e/p + Kp 2 )\ k (M) + C n p\ k {Mf' 2 . 



Therefore 



|A fc (r)-A fc (M)| < CmAs/p + p) 



provided that p < C 



As a corollary, for every fixed k we have Afe(T) — > Afc(M) as p — > and ^ — > 
and the convergence is uniform over all M E Ad. 
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The estimate (|1.4|) is a combination of Proposition l4~4l and Proposition l6. 31 where 
we prove an upper and a lower bound, respectively, for A&(F) in terms of Afc(Af). 
These propositions also provide somewhat sharper estimates on Afc(r) — Afc(-M). 
The second assertion of Theorem [T] follows from (|1.4[) and the fact that for every 
fixed k the eigenvalue Afc(M) is uniformly bounded over M G M. 

Our next result establishes convergence of eigenfunctions. Namely, it is possible 
to approximate an eigenfunction of Am corresponding to an eigenvalue A by a linear 
combination of eigenfunctions of Ar corresponding to eigenvalues close to A. The 
precise formulations are given in Theorems [3] and HI in Section [7J Here we give only 
a special case of this result where A has multiplicity 1. 

Theorem 2. Let f k be a unit-norm eigenfunction of —A m corresponding to an 
eigenvalue \ k — Afc(M) of multiplicity 1 , and let 5\ — min{l, A^+i — Afc, X k — Afc_i}. 
Then, for sufficiently small p and e/p (more precisely, if p + e/p < Cj^ k 5\), the 
eigenvalue Afc(L) of — Ap also has multiplicity 1, and for a corresponding unit-norm 
eigenvector u k we have 

\\Pfk-u k \\ L 2 (x) < C M ,k$x 1 (£/p + p), 

\\iuk - fk\\mM) < c M ,ko~x 1 ( £ /p + p)> 

where the norm \\ ■ \\l 2 (x) * s defined in ()1.2|) and the maps P: L 2 (M) — > L 2 (X) and 
I: L 2 (X) — > C 0,1 (M) are discretization and interpolation defined in Definitions ^. 1\ 
and \6.1[ respectively. 

Theorem [2] is a special case of Theorem |4l which handles arbitrary multiplicity. 
Theorem [3] is another variant where an estimate is uniform over M. (in particular, 
it does not depend on the size of spectral gaps) . However the rate of convergence 
guaranteed by Theorem [3] is worse than that in Theorem 

Remarks on the proof. Let us note that the upper bound 
(1.5) limsupA fc (r) < A fc (M) 

is nearly trivial. It follows from the fact that our graph Laplacian approximates the 
function Laplacian for every smooth function. Indeed, let /i, . . . , fk be orthonormal 
eigenfunctions of —Am with eigenvalues Ai(M), . . . , Afc(M). It is well known that 
the eigenfunctions are smooth (more precisely, their Cj* norms are bounded by 
CM,k, see [2])- Let u\, . . . , u k S L 2 {X) be discretizations of /i, . . . , f k . (For smooth 
functions the precise definition of discretization does not really matter; one can 
define e.g. Uj(x\) = fj(xi).) Since the functions fj are smooth, their discrete 
functions Uj almost orthonormal in L 2 (X) and their discrete Laplacians Ar(uj) 
are pointwise close to the Laplacians A^f/j - Hence 

(-A r Uj,Uj) w (-A M fj,fj)L3 = Xj(M) 

and therefore (— Aru,w) < Afc(M) for every u from the linear span of u%, . . . ,u k - 
Thus we have a fc-dimensional subspace of L 2 {X) where the norm of the discrete 
Diriclet energy functional (I2.4[) is bounded by approximately Afc(-M). By the mini- 
max principle it follows that Afc(F) < Afe(M), in other words, (|1.5[) holds. 

The proof of the upper bound in Sections |3j and [4] is different. We define a 
discretization map P: L 2 (M) -> L 2 (X) that makes sense for non-smooth func- 
tions and show that this map almost preserves the L 2 norm and almost does not 
increase the Diriclet energy, on a bounded energy level (see Definition 14.11 and 
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Lemma 14. 3p . This argument does not require pointwise eigenfunction estimates 
and yields sharper inequalities. 

The lower bound (i.e., the inequality liminf A/c(r) > Afc(M)) is more delicate. 
Here good approximation of Laplacians of smooth functions is not sufficient. For 
example, consider a disjoint union T of two graphs Ti and T2 each of which pro- 
vides a good approximation of the function Laplacian. The graph Laplacian Ar 
approximates the function Laplacian as well as Ar 1 and Ar 2 do, but the spectrum 
is different: every eigenvalue appears twice. 

To prove the upper bound, we construct a map /: L 2 (X) — > C 0,1 (M), called the 
interpolation map, with properties similar to those of the discretization map P, see 
Definition 16.11 and Lcmma l6.2l This map is essentially a convolution with a certain 
kernel (the form of the kernel is essential for the estimate in Lemma HT27 2)). With 
this map, the proof of the lower bound is similar to that of the upper bound. 

In addition, the maps P and / are almost inverse to each other on bounded 
energy levels (Lemma I6.4|) . These properties of P and / imply our eigenfunction 
estimates (Theorems [21 [3] and SJ by means of linear algebra arguments. 

Remark 1.3. The input data to the construction are p > and the finite metric 
measure space (X, p) . One naturally asks how sensitive are the resulting eigenvalues 
Afe(r) to "measurement error" in these data. A small relative error in weights pi 
results in a relative error of the same order in the L 2 (X) and the discrete Dirichlet 
energy (|2.4[) and hence to the eigenvalues. A small (of order e) variation of distances 
in X changes the set of edges of T: some edges of lengths p ± e may be added 
or removed. The discrete Dirichlet energy and hence the eigenvalues are clearly 
monotone with respect to adding edges. Therefore the eigenvalues are bounded 
above by those of the proximity graph defined by the parameter p 4- e in place of 
p, up to a factor (1 + e/p) n+2 . A similar argument yields a lower bound. 

Remark 1.4. The convergence of eigenvalues in Theorem [T] is uniform on a larger 
class of n-manifolds, namely those with bounded Ricci curvature and diameter and 
injectivity radius separated from zero. Indeed, by [1] this class is pre-compact in 
C 1 '" (and hence Lipschitz) topology. This pre- compactness and convergence for 
every individual manifold implies uniform convergence on the class. This can be 
shown with an argument similar to one outlined in Remark 11.31 

Remark 1.5. If the weights pi are constant (i.e., pi = po := vo\(M)/N), then the 
edge weights in our construction are also constant. Hence the graph Laplacian given 
by (|1.3|) is the ordinary (unweighted) graph Laplacian multiplied by a constant. 
Also note that in this case the degree in the graph is almost constant (up to a small 
relative error): the degree of every vertex approximately equals v n p n / Vo- 

Unweighted graph Laplacians has been studied much more thoroughly than 
weighted ones. If necessary, one can make the weights constant (at the expense 
of increasing the number of vertices) as follows. First approximate the weights pi 
by rational multiples of vol(M) and let q be a common denominator of these ra- 
tionals. Then replace every point Xi with weight pi = y vol(M) by pi points (at 
almost the same location) with weights equal to vo\(M)/q. The resulting metric 
measure space approximates (M, vol) as well as the original one do. 

Remark 1.6. Although our point is to avoid triangulation of a manifold, let us 
mention that triangulation-based techniques allow one to handle differential form 
Laplacians as well, see [8]. It is interesting whether a suitable generalization of a 
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graph Laplacian can be used for this purpose too. One can show that the spectrum 
of the differential form Laplacian is continuous with respect to Gromov-Hausdorff 
topology on M.. Hence a Gromov-Hausdorff approximation of a manifold (such 
as an e-net) determines differential form Laplacian eigenvalues up to a small error. 
However an explicit procedure of such determination is yet to be found. 

Organization of the paper. In Section [2] we collect various preliminaries. In 
Section [3] we prove some technical results about average dispersion in r-balls of 
a function / £ L 2 (M). This quantity, denoted by E r (f), is used throughout the 
paper as an intermediate step between Dirichlet energy in H 1 (M) and its discretiza- 
tion. In Section [4] we define the discretization map P and prove an upper bound 
for the graph eigenvalues (Proposition I4.4j) . Section [5] is devoted to properties of a 
smoothcning operator (the convolution with a special kernel) used in the definition 
of the interpolation map /. The key result there is Lemma 15.51 In Section [5] we 
define / and prove a lower bound for the graph eigenvalues (Proposition 16.31) . The 
proof of Theorems is contained in Section[7] In Section[8] (which is formally indepen- 
dent of the rest of the paper) we discuss various aspects of volume approximation 
in the sense of Definition 11.11 
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2.1. Discrete differential. Let E = E(T) be the set of directed edges of our 
graph. (Each pair of adjacent vertices gives rise to two elements of E.) Recall that 
every edge e,j = (xi,Xj) is equipped with a weight w(eij) = Wij. By L 2 (E) we 
denote the space of real- valued functions on E equipped with the following inner 
product: 



For a discrete function u : I-}Iwe define its discrete differential Su: E —¥ M. by 



in particular, (Aru,v) = (u, A-pv) and (ArU,u) — — \\Su\\ 2 for all u,v £ L 2 {X). 
(Here and almost everywhere in the paper we omit indices L 2 (X) and L 2 (E).) 
Thus Ar is self-adjoint and nonpositive on L 2 (X). 



2. Preliminaries 




(2.1) (6u)(eij) = u{xj) - u{xi). 

The discrete Dirichlet energy functional of T is the quadratic form 



(2.2) ui-» \\5u\\ 2 L2{E) = (Su,Su) L 2 (E) 
on L 2 {X). A straightforward calculation shows that 

(2.3) (A r u,-y) L 2 (x) = -(5u,8v) L 2(e), 
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The above consideration does not depend on a particular choice of weights. In 
our case, the discrete Dirichlet energy \\Su\\ = \\Su\\i,2tE) is given by 

(2.4) \\Su\\ 2 = — wKi,)-^)| 2 - 

P 

Since the operator — Ar is self-adjoint on L 2 (X) and the associated quadratic 
form is the discrete Dirichlet energy, the minimax principle applies: 

. (r . . \M\h(E) 

\k (I ) = mm max — — 

l ueL\o \\u\\ z L2{x) 

where the minimum is taken over all /c-dimensional subspaces L C L 2 {X). 

2.2. Local Riemannian geometry. Throughout the paper, M is a compact Rie- 
mannian manifold (without boundary) and n = dim M . The absolute values of 
sectional curvatures of M are bounded above by K and the injectivity radius is 
bounded below by Our standing assumptions are that p < io/2 and Kp 2 < 1/n 2 . 

For x £ M, exp x : T X M — > M is the Riemannian exponential map. We always 
restrict exp x to the ball £?2 P (0) C T X M, this restriction is a diffcomorphism onto the 
geodesic ball B 2p (x) and hence its inverse exp - ; 1 : B 2p {x) — > B 2p (0) is well-defined. 
We denote the Jacobian of exp^ at v € i?2p(0) C T X M by J x {v). 

By the Rauch Comparison Theorem, the relative distortion of metric by exp^, at 
v € B 2p (0) C T X M is bounded by 0(K\v\ 2 ) and hence 

(2.5) (1 + CnKlv] 2 )- 1 < J x {v) < 1 + CnK\v\ 2 . 
It follows that vol(i3 r (a;)) ~ v n r n as r — > 0, more precisely, 

|vol(S r (a;)) - v n r n \ < CnKr n+2 

for all r < 2p. 

The inner product in T X M defined by the Riemannian structure is denoted by (, ) . 
This scalar product allows one to identify T X M and T*M and we sometimes assume 
this identification to simplify notation. By grad/(x) we denote the Riemannian 
gradient of a function /: M — > M at x € M, i.e., the vector in T X M corresponding 
to the differential d x f 6 T*M . Recall that the gradient of the distance function 
d(-,y) at x is the velocity vector at the endpoint of the minimal geodesic from y to 
x, that is, 

(2.6) gradrf(-,y)( a; ) = - e ^ 1 y . 

2.3. Integration over balls. In this section we justify the normalization constant 
in (jl.3p . If Q is a quadratic form on R™ , then for every r > we have 

f v. r n+2 

(2.7) / Q(x)dx = — — -trace(Q). 
JB r (Q) « + 2 

Indeed, since both sides are preserved under orthogonal transformations and linear 
in Q, one can replace Q by its average under the action of the orthogonal group. 
Thus it suffices to verify (|2.7[) only for rotation-invariant quadratic forms, or, equiv- 
alently, for the form Q(x) = \x\ 2 . For this form, one computes the integral using 
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spherical coordinates: 

\x\ 2 dx = I t 2 Yoln^id B t {Q))dt = I nu n t n+1 dt- 



'Br(O) 

The identity ()2.7[) follows since trace(a; i-> \x\ 2 ) = n. 

Let / : R™ — > R be a smooth function. Integrating the Taylor expansion of / at 

x £ 1", 

f(x) - f(x ) = L(x - x ) + Q(x - x ) + o(\x - x Q \ 2 ), \x - x \ -t 0, 
where L = d Xo f and Q = ^d 2 f, using (|2.7|) . yields 

n+2 n-t-2 

(f(x)-f(x ))dx = ^_trace(Q)+ (r"+ 2 ) = ^—Af(x ) + o(r"+ 2 ) 
s r (^ ) n + 2 2(n + 2) 

as r -> 0. For a smooth function /: M — > R this relation holds as well since the 
Jacobian of the exponential map introduces an error term of order 0(r™ +3 ), as 
follows easily from (|2.5[) . Thus 

2{U t+2 I (f(x)-f(xo))dx^Af(x ) asp^O 

VnP n+2 J Bp (x ) 

for every smooth /: M — > R and every x G M. Furthermore the error term is 
controlled by the modulus of continuity of the second derivative of /. 

Replacing the above integral by the sum from (|1.3I) essentially replaces the inte- 
gration over the ball by integration over the union of the sets Vi (see Definition II. ip 
such that the respective points Xi belong to the ball. One easily sees that the error 
term introduced by this change is controlled by e/p 2 . (This estimate can be im- 
proved by introducing edge weights as in Remark ll.2j) . It follows that the discrete 
Laplacian of a smooth function approaches its ordinary Laplacian as p + e/p 2 — > 0. 

This observation is important for motivation of our definitions, but we do not 
use it in the proofs. Our arguments are based on the discrete Dirichlet energy and 
the minimax principle, these tools provide better estimates. 

3. Some estimates 

In this section we prove some inequalities for functions on M not involving 
discretization. 

Definition 3.1. Let / £ L 2 (M) and < r < 2p. For every measurable set V C M, 
define E r (f,V) £ R+ by 

E r (f,V)= [ [ \f{y)~f{x)\ 2 dydx. 

JV JB r (x) 

Let E r (f)=E r (f,M). 

Remark 3.2. The quantity E r (f) is bounded in terms of ||/||l2, namely 

(3.1) E r (f) < Cu n r n \\ff L2 . 

Indeed, 



E r (f)<2[ [ (\f(x)\ 2 + \f(y)\ 2 )dydx 

JM JB r (x) 

= 4/ / \f(x)\ 2 dydx = 4 f vo\(B r (x))\f(x)\ 2 dx, 

JM JBt(x) JM 
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and the right-hand side is bounded above by Cv n r n || 2 • Since E r is a nonnegative 
quadratic form, (|3.ip implies that it is a continuous map from L 2 {M) to R + . 

Lemma 3.3. Let f e H X {M) and < r < 2p. Then 

E r (f)<(l + CnKr 2 )^-r"+ 2 \\df\\ 2 L2 . 

n + 2 

Remark. The inequality turns to almost equality if / is smooth and r is small. 
This follows from the fact that the constant is the integral of the square of a 
coordinate function over the unit ball in R™, see (|2.7[) . 

Proof of Lemma HOI Since smooth functions are dense in H 1 (M) and E r is a con- 
tinuous map from H 1 (M) to R+, we may assume that / is smooth. Thus we can 
speak about pointwise values and derivatives of /. 
For every x G M , we have 

\f(y)-f(x)\ 2 dy= [ \f(exp x (v))-f(x)\ 2 J x (v)dv 

>B r (x) JB r (0)CT x M 

where J x is the Jacobian of exp^, (see Section [2. 2 p and J dv denotes the integration 
with respect to the Euclidean volume on T X M determined by the Riemannian scalar 
product. Since J x (v) < 1 + CnKr 2 for all v £ B r (0) C T X M, it suffices to prove 
that 

(3.2) A:= [ I |/(exp x («))-/(x)| 2 ^<-^-r"+ 2 ||d/||i 2 . 

J M J B r {0)<ZT x M n + Z 

For every x and v we have 

/(exp»)-/0r)= f £j(en> x (tv))dt= [ #($ t (x,v))dt 
Jo Jo 

where <I> t : TM —> TM is the time t geodesic flow, namely $*(x, v) = {jx,v{t),j' x v (t)) 
where j XjV is the constant-speed geodesic given by j x>v (t) = exp x (t«). In the ex- 
pression df(<& t {x, v)), the derivative df is regarded as a (fiberwise linear) map from 
TM to R. 

The above identity and the Cauchy-Schwartz inequality imply that 

|/(exp») - f(x)\ 2 < f \df(1> t (x,v))\ 2 dt. 
Jo 

Hence the right-hand side of (|3.2p can be estimated as follows: 

(3.3) A< [ f \df(^ t (0)\ 2 dvol T M(Odt 

Jo JB(r) 

where B(r) C TM is the set of all tangent vectors £ £ TM such that |£| < r, and 
voItm is the standard 2n-dimensional volume form on TM. Since B(r) is invariant 
under $ t and $t preserves voItm (by Liouville's Theorem), the inner integral in 
p.3|) does not depend on t. Therefore 

A< f \df(0\ 2 dvo\ TM (0 = [ I \d x f(v)\ 2 dvdx 

JB(r) Jm J B r (0)cT x M 



where the second identity follows from (|2.7j) . This proves (J3T2J) and hence the 
lemma. □ 
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Lemma 3.4. Let < r < 2p, f € L 2 (M) and V C M be a measurable set such 
that vol(V) = n > and diam(V") < 2e where e < r. Let a = yT x J v fix) dx be the 
integral mean of f\v ■ Then 

j v \m-afdx<-^E r{ f,V). 
Proof. A standard computation shows that 

(3.4) J | /(*) ~ a \ 2 dx=±-J^ J \f(x) - f(y)\ 2 dxdy. 

Fix x, y € V and consider the set U = B r (x) n B r (y). Observe that U contains the 
ball of radius r — \xy\/2 > r — e centered at the midpoint between x and y. Hence 
vo\(U) > Cv n (r — e) n . For every z g U we have 

\f(x) - f(y)\ 2 < 2(1/(1) - f(z)\ 2 + |/(2/) - /(z)| 2 ). 

Therefore 

\f(x) - f{y)\ 2 < Jj\f{x) - f(z)\ 2 + \f(y) - f(z)\ 2 ) dz 

< + ^)) < (FW + ^)) 

where 

F(x)= / |/(x)-/(z)| 2 dz. 

iB r (i) 

Plugging the last inequality into ()3.4[) yields 

X I'M - »' 2 * s ^F^F X X {F{X) + F{V)) dxi » 

The lemma follows. □ 

4. Discretization map and upper bound for Afc(r) 

Let X — {xi}fL 1 C M and /i be as in Theorem [TJ Recall that p is a measure 
on X and (X, fi) e-approximates (M, vol) in the sense of Definition 11.11 We fix 
a partition {Vi}fL 1 of M realizing this approximation, that is, Vi C B e (xi) and 
vol(Vi) = fii := n(xi) for each i. We assume that e < p/n. 

Definition 4.1. Define a discretization map P: L 2 (M) — > L 2 (X) by 

Pf(x l )= f i- 1 f f(x)dx. 

JVi 

In other words, Pf(xi) is the integral mean of /|v;. 
We also need a map P* : L 2 (X) -> L 2 (M) defined by 

TV 

P*u = $^u(a;<)lv; 
i=i 

where 1\a is the characteristic function of the set V,. 
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From the Cauchy-Schwartz inequality one easily sees that 
(4-1) \\Pf\\ < ll/IU* 

for every / £ L 2 (M), where the norm in the left-hand side is defined by (|1.2|) . The 
definition implies that P* preserves the norm: 

\\P*U\\ L * = \\U\\ 

for all u £ L 2 (X), and is adjoint to P: 

{f,P*U) L 2( M ) = (Pf,u) L 2 {X) 

for all u £ L 2 (X), f £ L 2 (M). 

Lemma 4.2. Let f £ H 1 (M). Then \\f - P*Pf\\ L 2 < Cne\\df\\ L 2. 
Proof. We have 

\\.f-P*Pf\\ 2 L 2=J2 ! \f{x)-Pf{x t )\ 2 dx. 

i Jv i 

By Lemma 13.41 for every r £ (e, 2p) and every i we have 

\f(x) - Pf{ Xl )\ 2 dx < ° E r (f, Vi). 
Vi v n [r — e) 

Note that £\ E r(f, Vi) = E r {f) by definition. Therefore 

C C r n 

||/ - P*Pf\\h < , , E r (f) < — — - r 2 \\df\\ 2 L - 2 



where the last inequality follows from Lemma 13.31 Now let r = (n + l)e, then 
(7 r 7F = ( 1 + ^)"< 3 ' hence 

\\f~P*Pf\\h < ^r 2 \\df\\ 2 L2 <Cns 2 \\df\\ 2 L2 . 
The lemma follows. □ 
Lemma 4.3. Let f £ H 1 (M). Then 

(1) |P7||-||/|U»| <Cne\\df\\ L 2; 

(2) ||<y(P/)|| < (1 + o-)\\df\\ L 2 where a = Cn{Kp 2 + e/p). 
Proof. (1) Since P* preserves the norm, we have 

IIIP/II - II/IM = \\\P*Pf\\L* - II/IU < II/- P*Pf\\ L * < Cne\\df\\ L 2 



where the last inequality follows from Lemma 14.21 
(2) By 424} we have 

V-n P 

The definition of Pf implies that 

Pf( Xj ) - Pfixi) = — [ [ (f(y) - f(x)) dydx. 

piPj JVi JVj 

Hence, by the Cauchy-Schwartz inequality, 

\Pf( Xj )-Pf( Xi )\ 2 <— f [ \f(y)-f(x)\ 2 dydx. 
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Therefore 



W)II 2 <^E E / / \f(y)-f(x)\ 2 dyd x 



n + 2 



v nP JM Ju(x) 



Vi J Vj 

\f(y)-f(x)\ 2 dydx 



where the set U{x) C M is defined as follows: if x € Vj, then U(x) = []j. x .^ x . Vj. 
Note that U(x) C B p+ 2 e (x). Hence 

77-1-9 

\\5{PfW<— ^E p+2s (f) 



E p+2 e{f) < ^(p + 2eT +2 (l + a 1 )\\df\\l 2 . 



By Lemma 13.3 



where o\ = CnKp 2 . Therefore 

||<5(P/)|| 2 < (1 + 2e/p) n+2 {l + CTl )||d/||| 2 < (1 + a)\\df\\l 2 . 
where a — Cn(Kp 2 + e/p). This finishes the proof of Lemma T4.3I □ 
Proposition 4.4. Let \ k = \ k {M), keN. Then 

A fe (r) < (l + S{e,p,X k ))X k 

where 

S(e, p, A) = Cn(Kp 2 + e/p + eVX), 
provided that £-\/Afc < c/n. Here C and c are absolute constants. 

Proof. By the minimax principle, it suffices to show that there exists a linear sub- 
space L C L 2 {X) such that dimL = k and 

\\Su\\ 2 

sup < (1 + S(s,p, Afc))Afe. 

ueL\{o} Nr 

Denote A = A fe . Let W c ^(M) be the linear span of orthonormal eigenfunc- 
tions of —Am corresponding to eigenvalues Ai, . . . , Afe. For every / £ VF, we have 
||d/||£ 2 < A||/|| 2 2 . By Lemma|13i;i) it follows that 

\\Pf\\ > ll/IU' - Cne||4f|U« > (1 - CneVA)||/|| La 

for every f € W. Hence P\w is injective if eV\ < l/Cn. Let L = P(W), then 
dimi = k. Pick u e L \ {0} and let / G W be such that u = Pf. Then 

and, by Lemma 



|M| 2 >(l-eWA)||/|| 2 2 



||fe|| 2 <(l + CTl )||d/||l 2 <(l + a 1 )A||/||l 2 
where o\ = Cn(Kp 2 + e/p). Hence 

||<5m|| 2 „ (1 + (Ti)A „ ,„ ... , 

and the proposition follows. □ 
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5. SMOOTHENING OPERATOR 

In this section we prepare technical tools for the interpolation map, which is 
defined in the next section. These tools are independent of the discretization. 
Define a function ip : R + — > R + by 

'^(i-* 2 ). 0<t<l, 
0, t>l. 



The normalization constant is chosen so that 



(5.1) / ip{\x\)dx = l. 
Indeed, by (|2.7|) we have 

ft I |2\ j nV " 

(1 — J dx = v n 



3,(0) n + 2 n + 2 

Fix a positive r < 2p and consider a kernel k r : M x M —> R + defined by 

k r (x,y) = r~ n ip(r~ 1 \xy\) 
and the associated integral operator A° : L 2 (M) — > C 0,1 (M) given by 

A° r f(x)= f f(y)k r (x,y)dy. 

J M 

Note that k r (x,y) — k r (y,x) and 

i. / \i Cn 
\K{x,y)\ < 



v n r- 

for all x, y £ M . A direct computation (using the derivative of the distance function, 
see (|2.6p ) yields 

n + 2 

(5-2) gr&dk r (-,y)(x) = — -exp" 1 ^) 

v n r + 

for y € B r (x). 

Define £ C° 4 (M) by 9 = A°(1 M ). If the metric of M were flat, we would have 
6 = 1m by (|5.ip . The following lemma estimates \\0 — 1m\\ m t ne Riemannian case. 

Lemma 5.1. For every x € M, one has 

(1 + CnATr 2 )" 1 < 6{x) < 1 + CnKr 2 

and 

\d x 9\ < Cn 2 Kr. 

Proof. By definition, 

8(x)=r~ n [ ^(r- 1 d(x,y))dy = r- n [ ijj{r- l \v\)J x {v) dv 

JB T (x) JB r (Q)C.T x M 

where J x {v) is the Jacobian of the Riemannian exponential map, see Section 12.21 
Since the integral of ip equals 1, the Jacobian estimate (|2.5p implies the first asser- 
tion of the lemma. 

To estimate \d x 9(x)\, we compute it using (|5.2[) : 

n + 2 f n + 2 f 
grad0(z) = — =• / exp~ 1 (y)dy= — j / vJ x (v)dv. 

V n T J B r (x) v nf JB r (0)CT x M 
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Since J B vdv — 0, one can replace J x {v) in the last integral by J x (v) — 1. Then 
the Jacobian estimate (12.51) implies that 

77 + 2 f 

\d x 6{x)\ = |grad0(a;)| < — -—^ / |v| • CnKr 2 dv < Cn 2 Kr 

(the last inequality follows from the relations |tj| < r and vol(i? r (0)) — v n r n ). □ 
Definition 5.2. Now we define a bounded operator A r : L 2 (M) — > C 0,1 (M) by 

A r / = e- 1 ■ A° r f. 

The factor # _1 ensures that A r preserves the subspace of constants. 
Lemma 5.3. For every f £ L 2 (M) one has 

||A r /|| L 2 < (l + <r)||/|| L2 

where a — CnKr 2 . 

Proof. This is a standard estimate. For every x € M we have 

2 



|A°/(z)| 2 = 



f(y)k r (x,y)dy 



M 



< [ l^k r (x,y)dy\ ^JJf(y)\ 2 k r (x,y)dy 



= 0{x) / \f(y)\ 2 k r (x,y)dy 

J M 

by the Cauchy-Schwartz inequality. Hence 

\A r f(x)\ 2 = 9(x)- 2 \A°J(x)\ 2 < Oix)- 1 f \f(y)\ 2 k r (x, y) dy 

< (1 + a) [ \f(y)\ 2 k r (x,y)dy 
by Lemma 15. II Integrating this inequality over M yields 



\\A r f(x)\\i 2 <(l + a) \f(y)\ z kr(x,v)dxdy<(l+*y\\f\\b 
Jm Jm 

since for every y € M we have 

k r (x,y) dx = 9{y) < 1 + a 



I M 

by Lemma 15. II The lemma follows. □ 
Lemma 5.4. For every f £ L 2 (M) one has 

(5-3) \\A r f-f\\ 2 L2 <-^E r (f). 

Proof. We fix a particular function f:M — > M. representing the given element of 
L 2 (M), so we can speak about pointwise values of /. Fix x £ M and let a = f{x). 
Since A r preserves the constants, we have 

Ar/(aO - f(x) = A r f(x) -a = A r (f - a ■ l M )(x) 
= 0-\x) f (f(y)~a)k r (x,y)dy 

J B T (x) 

= 9-\x) f (f(y)-f(x))k r (x,y)dy. 

J B r {x) 
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By the Cauchy-Schwartz inequality, it follows that 

\A r f( X )-f(x)\ 2 <0- 2 (x)([ k r (x,y)dy]([ \f(y)-f(x)\ 2 k r (x,y)dy) 

\JB r (x) J \JB r {x) J 

= 0- 1 {x)[ \f(y)- f(x)\ 2 k r (x,y)dy 

J B r (x) 

<-^f \f(v)-f(-)\ 2 

v n r n J Br ( x ) 

since |6 l_1 (x)| < C (cf. Lemma |5~T|) and \k r (x,y)\ < "„ . Integrating this inequal- 
ity yields the result. □ 

Lemma 5.5. For every f G L 2 (X) one has 

n 4-9 

\\d{Kf)\\h < (l + Cn 2 Kr 2 )—^E r (f). 

Proof. We fix a particular function / : M — > R representing the given element of 
L 2 (M), so we can speak about pointwise values of /. Denote A r f by /. For any 
constant a € R we have 

f(x) =a + e-\x) f (f(y)-a)k r (x,y)dy 

JB r (x) 

for every x € M. Differentiating this identity yields 
dj=d- 1 (x)[ (f(y)-a)d x k r (;y)dy + d x (8~ 1 ) f (f(y)-a)k r (x,y)dy. 

JB r {x) JB T {x) 

Substituting a = f(x) yields 

d x f ^e- 1 (x)A 1 (x)+A 2 (x) 

where 

Ai(x)= f (/(»)- f(x))d x k r (-,y)dy 

JB r (x) 

and 

A 2 (x) = d^Q- 1 ) f (f(y) ~ f(x))k r (x, y) dy. 

JB r (x) 

Since |(9 -1 (a;)| < 1 + CnKr 2 for all x £ M (cf. Lemma f5 . 1 [) . we have 
(5.4) \\df\\ L 2 = + A 2 \\ L 2 < (1 + CnKr^WAxWv + \\A 2 \\ L 2. 

Let us first estimate || ^4.2 Hx- 2 - By the Cauchy-Schwartz inequality, 

\A 2 (x)\ 2 < {d^O-^l 2 I f k r (x,y)dy)lf \f(y)-f(x)\ 2 k r (x,y)dy) 

\JbAx) J \JB r (x) J 

= \d x (6- 1 )\ 2 9(x) [ \f(y)~f(x)\ 2 k r (x,y)dy 

J B r {x) 

\f(y)-m\ 2 dy 

v n r J B T {X) 
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since 9(x) < C, ^(fl" 1 )! < Cn 2 Kr (cf. LemmaP and \k r (x,y)\ < Inte- 
grating this inequality yields 

We rewrite this inequality as follows: 



(5-5) \\A 2 \\ L2 <Cn 2 Kr 2 ^^-^E r {f). 

Now let us estimate A\. Fix x € M. Recall that 

\Ax{x)\ = vaax{{Ai{x),w) : to £ T X M, |w| = 1} 

where the angle brackets denote the standard pairing of co- vectors and vectors. Let 
to G T X M be a unit vector realizing this maximum. Then |Ai(x)| = (Ai(x),w). 

Plugging the expression (15.21) for the derivative of k r into the definition of Ai 
yields 

n + 2 f 

\A 1 (x)\={A 1 (x),w) = / (f(y)- f(x))(ex Px 1 (y),w)dy 



n + 2 



B r (x) 

tp(v)(v,w)J x (v)dv 



V n r n+2 J B r (0)CT,M 

where ^j(u) = f(exp x (v)) — /(a;). For brevity, we denote the ball B r (0) c T X M 
by B. By the Cauchy-Schwartz inequality, it follows that 

n + 2 x " 
v n r 

Since |w| = 1, we have 

m -i- 9 r 

(v, w) dv = 1 



|^i(^)| 2 < ( 7—^+2 ) ( / WtUvfdv) ( / (v,w) 2 dv 



v n r n+2 

hence the above inequality boils down to 



n + 2 
n + 2 



\Arix)? < -f-^ / \ V (v)fj x (vfdv 



n+2 



\f(y)~f(x)\ 2 J x (exp-\y))dy 

B r (x) 

<{ l + CnKr 2 )^- 2 I \f(y)-f(x)\ 2 dy 
v n r JB r (x) 

where the last inequality follows from the Jacobian estimate (|2.5|) . 
Integrating this inequality with respect to x over M yields 

\\Mh<(l + CnKr*)-^E r (f). 
This, and (I53|) imply that 

\\df\\v < ((1 + CnKr 2 f 2 + Cn 2 Kr 2 ) J J^±E r (f) 

<(i + ^ 2 ^ 2 )/|^C0- 

The lemma follows. □ 
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6. Interpolation map and lower bound for Afe(r) 
Definition 6.1. Define the interpolation map I: L 2 (X) — > C 0,1 (M) by 

Iu = A p - 2e P*u 

where A p _2 £ is the smoothening operator defined in the previous section (see Defi- 
nition EH]) and P* : L 2 (X) -> L 2 (M) is defined in Definition 



Lemma 15.31 and the fact that P* preserves the norm imply that 
(6.1) ||iu|U2 < (1 + CnKp 2 )\\u\\ < C\\u\\ 

for every u € L 2 {X). 

Lemma 6.2. For every u 6 L 2 (X) one has 

(1) |||/«|| La -||«|||<Cp||HI; 

(2) \\d(Iu)\\ L 2 < (1 + <t)\\Su\\ where a = Cn 2 Kp 2 + Cns/p. 
Proof. (1) Since ||P*u||l2 = ||u||, we have 

|||/u|| ia - ||u||| = \\\Iu\\ L 2 - \\P*U\\ L 2\ < \\Iu-P*u\\ L 2 



By Lemma [5T41 

\\Iu-P*u\\ 2 L2 = \\A p ^ E P*u-P*u\\ 2 L2 < —^—- E p . 2s (P*u) 
Since e < p/n and e < p/3, we have (p — 2e)~" < Cp~ n , hence 
(6.2) ||Ju - P*u\\ 2 L2 < ^- Ep ^ £ (P*u). 

Let us estimate P p _2 £ (P*'u) in terms of 5u. By definition, 

2 _ n + 2 



HI 2 = ^j\u(Xj) - u(Xi)\ 2 

Vfl p 

n+ ^2 [ I \P*u{y)-P*u{x)\ 2 dydx 
jm Ju(x) 



where sets U(x) C M are defined as follows: if x g Vt, then f/(x) = {Jj. x .^ x . Vj. 
Since E/(x) D P p _2 £ (:e), we have 

||H| 2 >^2-/ / \P*u{y)-P*u{x)\ 2 dydx = ^±^E p . 2e {P*u). 

VnP n+2 Jm Jb o - 2c (x) VnP n+2 



Thus 

n + 2 



n+2 
n " II X II 2 



(6.3) P p -2e(P*^) < -rr^-ll^ 



This and (jOf imply that 

z/ n p™ n + 2 



Iu P*u\\ 2 L2 < ^-^-\\Su\\ 2 < Cp 2 \\6u\\ 2 . 



Hence 

(6.4) \\Iu-P*u\\ L 2 < Cp\\6u\\ 

and the first assertion of the lemma follows. 
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(2) By Lemma E3 

77+2 

\\d(Iu)\\% = \\d{k p _ 2E P*u)\\l 2 < (l + ai ) —^—L-^ E p - 2e (P*u) 
where a\ = C'n 2 Kp 2 . By 

77-1-9 / \ n 

Un{p X)^ Ep -^ {p * u) -h^) IIHI 2 <(i + ^)||H 2 - 

where a% = Cue/ p. Thus 

\\d(Iu)\\ 2 L3 < (1 + (Ji)(l + cr 2 ) ||^-«|| 2 < (1 + <t)\\Su\\ 2 
where a — Cn 2 Kp 2 + Cne/p. The second assertion of the lemma follows. □ 

Proposition 6.3. Let X k = X k (M), fceN. Then 

A fe (r) > {l-S(e,p,X k ))X k 

where 

6(e, p, X) = C(n 2 Kp 2 + ne/p + p\f~X) 
provided that p\fXk < cq . Here C and cq > are absolute constants. 

Proof. With Lemma 16.21 at hand, the proof is similar to that of Proposition 14.41 
Let A = Afc(r). We assume that A < X k (M), otherwise there is nothing to prove. 
By the minimax principle, it suffices to show that there exists a linear subspace 
L C H 1 (M) such that dimL = k and 



sup 



WW 



L 2 



< (l-J(e,p,A)) _1 A. 



feL\{o} 11/11 

Indeed, this inequality would imply that X k (M) < (1 — 5(e, p, A)) _1 A and therefore 

A > (1 - *(e, p, X))X k > (1 - S(e, p, X k ))X k . 

Let W C L 2 (X) be the linear span of k orthonormal eigenvectors of — Ar cor- 
responding to eigenvalues Ai(r), . . . , Afc(F). For every u € W, we have ||5it|| 2 < 
A||w|| 2 . By Lemma IQl) it follows that 

\\Iu\\ L 2 > ||u|| - Cp||«5«|| > (1 - C P VX)\\u\\ 

for every u <E W. Hence I\w in injective if CpvX < 1. Let L = I(W), then 
dimi = k. Pick f e L \ {0} and let u £ W be such that / = Iu. Then 

\\ff L2 >(i-CpVx)\\uf 

and, by Lemma E2JJ2), 

\\df\\l 2 <(l + a 1 )\\Su\\ 2 <(l + a 1 )X\\u\\ 2 . 
where a\ = 1 + Cn 2 Kp 2 + Cne/p. Hence 



\\df\\ 2 L2 (1 + qQA 
||/|| 2 L2 ~ 1-CpVX 



» <^±^< {1 - S (e,p,X))-iX 



and the proposition follows. □ 

We conclude this section by showing that the operators P and I are almost 
inverse to each other (at bounded energy levels). 
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Lemma 6.4. 1. Let f e H 1 (M). Then 

\\IPf-f\\ L 2<Cp\\df\\ L2 . 

2. Let u £ L 2 (X). Then 

\\PIu-u\\ < Cp\\Su\\. 
Proof. 1. Let f = P*Pf. By definition, 

\\IPf - f\\l? = \\Ap-2ef ~ f\\v < l|Ap-2e(/ " f)\\» + Up-2ef ~ f\\&. 

Lemma 15.31 and Lemma 14.21 imply that 

||A,_a B (/- < CWJ-fWv < Cne\\df\\ L 2. 

Next, by Lemma T5. 41 and Lemma T3.31 

(In 

||A p _ 2s / - /||| 2 < Vn{p _ 2e)n Ep-2e{f) < Cp 2 \\df\\\ 

Combining the above inequalities and using the fact that e < p/n, we obtain the 
first assertion of the lemma. 

2. Since P* preserves the norm, we have 

\\PIu - u\\ = \\P*(PLu - u)\\ < \\P*PIu - Iu\\ L 2 + \\Iu - P*u\\ L 2. 

By Lemma 14.21 

\\P*PIu - Lu\\ L 2 < Cne\\d(Iu)\\ L 2 < Cne||<5u||, 
where at the last stage we use Lemma I6T2T 1). By (|6.4|) . 

\\Iu-P*u\\ L 2 < Cp\\8u\\ 
As e < p/n, the above inequalities imply the second assertion of the lemma. □ 

7. ElGENFUNCTION APPROXIMATION AND PROOF OF THEOREMS 

To prove Theorem HJ first observe that the estimate (|1.4p follows from Proposi- 
tions 133] and [231 Next recall that, as follows from [TJ, the space A4 is pre-compact 
in Lipschitz topology. Therefore the eigenvalue Afc (M) is uniformly bounded for all 
M € M., that is, A/-(M) < Cjw,k- Using this fact, we obtain the second estimate in 
Theorem Q] from the first one. 

To proceed with the eigcnfunctions approximations, we introduce some notation. 
For an interval Jcl, denote by Hj(M) the subspace of H 1 (M) spanned by the 
eigenfunctions with eigenvalues from J. In particular, -H{a} (M) is the eigenspace 
associated to an eigenvalue A. We abbreviate i?(_ oo A )(M) by H\{M). We use 
similar notation Hj(X) and H\(X) for subspaces of L 2 (X) spanned by eigenvectors 
of -A r . 

Note that the dimension of H\(M) is uniformly bounded over M G M (for every 
fixed A), see [51 Theorem 3]. 

We denote by Pj the orthogonal projector from L 2 (M) to Hj(M). Note that 
Pj does not increase the Dirichlet energy norm. Similarly to the above notation, 
we abbreviate P(-oo.a) by F\. We use the same notation Pj and Pa for orthogonal 
projectors from L 2 (X) to Hj(X) and H\{X). 
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Lemma 7.1. 1. Let A > and f e H\{M). Then 

\\5{Pf)\\>(l-a)\\df\\ L , 

where a = C{p\TX + n 2 Kp 2 + ne/p). 
2. Let A > and u e H\(X). Then 

\\d(Lu)\\ L2 > (l-a)\\6u\\ 

where a — C(p\/~X + nKp 2 + ne/p) 

Proof. 1. First we are going to estimate \\d(I P f)\\ ^2 from below in terms of ||d/||i2. 
Since the projector Pa : L 2 (M) — > H\(M) does not increase the Dirichlet energy, 

\\d(IPf)\\ L 2 > \\d(F X IPf)\\ L Z > \\df\\ L 2 - \\d(P X LPf-f)\\ L 2. 

Since / e H\(M), we have 

\\d(P x IPf - f)\\ L3 = \\d(P x (IP.f - f))\\ L2 < V\]\¥x(IPf - /)|U= 

< V\\\IPf-f\\ L 2<Cpy/X\\df\\ L 2. 

where the first inequality follows from the fact that ||d<7||j,2 < vA||5||x,2 for every 
g E H\[M) and the last one from Lemma \6. 41 Thus 

\\d(LPf)\\ L 2 > (l-aOHd/H^ 

where o\ = Cp\/\. By Lemma 

||<IP/)|| L2 <(l + a 2 )||5(P/)|| 

where er 2 = C(n 2 Kp 2 + ne/p). Thus 

\\S(Pf)\\ > (l + a^il-a^WdfW^ 

and the first assertion of the lemma follows. 

2. The proof of the second assertion is completely similar. Just interchange the 
roles of M and X and use Lemma I4.3f 2) rather than Lemma I6.2f 2) at the final 
step. □ 

We need the following simple lemma from linear algebra. 

Lemma 7.2. Let L be a finite- dimensional Euclidean space and k = dimL. Let Q 
and Q 1 be quadratic forms on L and Ai < • • • < Afc and X'i < . . . X' k their respective 
eigenvalues. Suppose that Q > Q' . Then 

sup {Q(v) - Q'(v)} < k max {Xj - A' }. 

«6i,||«||=l l<j<k 

Proof. The left-hand side is the largest eigenvalue of the quadratic form Q — Q' . 
Since Q — Q' is nonnegative, its largest eigenvalue is bounded above by it trace. 
On the other hand, 

k 

trace(Q — Q') = trace(Q) — trace(Q') = /^(Ao — A') < k max {A,- — A'}, 

hence the result. □ 

We fix orthonormal eigenfunctions {fk}kLi of —Am and orthonormal eigenvec- 
tors of -A x . 
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Lemma 7.3. 1. Let A = X k (M). Then for every a > 0, 

\\Pf k - Px+aPfkf < C MM a-\p + Elp) 

and 

\\S(Pfk -^x+aPfk)\\ 2 < C Mtk {l + a- l )(p + s/p) 
provided that p + e/p < C7} k . 

2. Let A = Xk(X). Then for every a > 0, 

||J-ujfe — PA+ a J«jfc||i= < C M<k cr l (p + s/p) 

and 

\\d{Lu k - V x +aIu k )\\ 2 L2 < CmM 1 + 0(P + e IP) 
provided that p + e/p < CJ^ k . 

Proof. 1. Let W be the linear span of fx, . . . , f k and L = P(W) C L 2 (X). As in 
the proof of Proposition 14. 4[ we have dim L = kiip + e/p is sufficiently small. Let 
Q denote the discrete Dirichlet energy form on L 2 (X), and let Af < ■ ■ ■ < X k be 
the eigenvalues of Q\l (with respect to the Euclidean structure on L defined by the 
restriction of the L 2 (X) norm). 

Recall that Afe(M) < Cm^. This and Lemma T4. 31 imply that for every / € W, 

(l-a)\\f\\ L2 <\\Pf\\ <(l + o-)\\f\\ L i 

and 

\\S(Pf)\\ <(1 +<T)\\df\\ L 2 

where a — CM.k{p + e l P)- By the minimax principle it follows that 

(7.1) \f < (i±2)\.(M) < \ 3 {M) + C M MP + e/p) 

for all j < k, provided that a < 1/2. 

Now define another quadratic form Q' on L 2 {X) by 

Q'{u) = Q(P x+a (u)) + X\\u- ¥x +a (u)\\ 2 . 

Clearly Q' < Q. The eigenvectors u\,ui,... of Q are also eigenvectors of Q' and 
the corresponding eigenvalues are Ai(T), A2(r), . . . , A m (r), A, A, ... , where m is the 
largest integer such that A m (T) < A + a. Therefore for every j < m and every 
j-dimensional subspace V C L 2 (X) we have 

(7.2) sup ^>min{A,A*(r)}. 

Indeed, V has a nontrivial intersection with the orthogonal complement of the 
linear span of u±, . . . , Uj—x, and any vector v from this intersection satisfies Q'[v) > 
min{A, Aj-(r)}||u|| 2 . Let A^ < • • • < A' fc be the eigenvalues of Q'\l (with respect to 
the restriction of the L 2 (X) norm to L). Then (|7.2[) and the minimax principle 
imply that X'j > min{A, Aj(T)} for all j < k. By Theorem [T] it follows that 

A;- > X j (M)-C M ,kip + e/p). 

and hence, by (|7.1|) . 



xf — x^ < CmAp + z/p) 

for all j < k. This and Lemma 17.21 imply that 

(7.3) Q{u) ~ Q'(u) < C M>k (p + e/p)\\u\\ 

for every u € L. 
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Let u £ L and vf = u — Pa+ q u- Since Q'(¥\+ a u) — Q(V\+ a u), we have 

(7.4) Q(u) - Q'(u) = Q(v!) - Q'{v!) = Q(u') - \\\u'\\ 2 > ^—Q(u') 

A + a 

where the last inequality follows from the fact that Q{u') > (A + a) || it' || 2 since 
v! e H[ X + a ,+oo){X). Now (fOl) and (|7Tl) imply that 

Q(u')<—C M . k (p + e/p)\W\\ 2 

a 

and therefore 

|| U '|| 2 < {X + a^Qiu') < a^CM^p + e/p^uf. 

Substituting u — P f k into the last two inequalities and taking into account that 
II P /fc|| < 1 + <t < 2 yields the first assertion of the lemma. 

2. The proof of the second assertion is similar. Just interchange the roles of M 
and X and use Lemma 16.21 rather than Lemma 14.31 □ 



Lemma 7.4. 1. Let A = Afe(Af) and let a, (3,^ > be such that a < (3 < 7 < 1 
and the interval (A + a, A + /?) does not contain eigenvalues of — Ar- Then 

\\Pf k - P (A _ 7iA+Q] P/ fe || 2 < Caj- 1 + CM^-^ip + e/p) 

provided that p + e/p < CT^ k . 

2. Let A = Afe(r) and let a, ,8,7 > be such that a < (3 < 7 < 1 and the interval 
(A + a, A + /?) does not contain eigenvalues of —Am- Then 

\\Iu k - P(\-y,x+ a ]Iu k \\ 2 L2 < Caj- 1 + CVfc/rV 1 ^ + f -Ip) 

provided that p + e/p < C7^ k . 

Proof. 1. As in the previous lemma, we denote the discrete Dirichlet energy form 
by Q. Let u = Pfk- Decompose u into the sum of three orthogonal vectors u — uq + 
+ where u G #(a- 7 ,a+ q ]PQ, U- € H^^x-j]^) and u+ G H {x+a , + oo)(X). 
Note that it + G Htx+p,+oo)(X) due to our assumption about eigenvalues of — Ar- 
Applying Lemma 17.31 with /3 in place of a yields that 

(7.5) \\u+\\ 2 ^Cm^-Hp + s/p) 

and 

Q(u+) < C M ^-\p + e/p). 

By Lemma I7TTI 

Q(u) = \\S(Pf k )\\ 2 > (1 - ^OIMAII^ - (1 - <7i)A 
where ci = C'M,k(p + z/p)- Therefore 

Q{u a ) + Q(u-) = Q{u) - Q{u+) > A — cr 2 
where 0-2 = CM,k/3^ 1 {p + s/p)- On the other hand, 

Q(u Q ) < (A + a)|| Mo || 2 



and 
hence 



Q(«-)<(A- 7 )IMI 2 , 

A - CJ2 < Q(«o) + Q(«-) < A(|M| 2 + \\u-\\ 2 ) + a\\u \\ 2 - j\\u-f 
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Observe that 

holl 2 < KII 2 + IMI 2 < IMI 2 = II^MI 2 <i + <T 3 

for 03 = CM.k(p + £/p), where the last inequality follows from Lemma EOT l). Thus 
A - o-2 < A(l + cr 3 ) + a(l + cr 3 ) - 7l|w-|| 2 , 

or, equivalent ly 

\\u-\\ 2 < 7~ V2 + Act 3 ) + a 7 - 1 (l + 0-3). 
The right-hand side is bounded by C M ,kl~ l fi~ l {p + e/p) + Coq~ x . This and (173)) 
yield the first assertion of the lemma. 

2. The proof of the second assertion is similar, with the roles of M and X 
interchanged. □ 

Theorem 3. 1. Let X = Xk(M) and let fk be a corresponding unit-norm eigen- 
function of —Am- Then for every 7 G (0, 1), 

\\Pfk -P (A - 7 ,A+ 7)J PM| 2 < C M ,kJ- 2 (p + £/p) 1/2 

provided that p + e/p < Cj^ k . 

2. Let A = Afe(r) and let Uk be a corresponding unit-norm eigenfunction of — Ap. 
Then for every 7 6 (0, 1), 

||/u fc -P (A _ 7iA+7) /« fc ||| a < C M ,kJ~ 2 (p + e/p) 1/2 

provided that p + e/p < Cj^ k . 

Proof. Plug a = ft = (p + e/p) 1 / 2 "/ into Lemma [7T41 Since the interval (A+a, A+/3) 
is empty, the assumption about eigenvalues is satisfied automatically. The desired 
estimates follows from Lemma [7.41 and the relations a < 7, wy^ 1 = (p + e/p) 1 / 2 
and/3- 1 = 7 - 1 (p + e /p)- 1 /2. □ 

The next theorem provides somewhat sharper estimates (which are however not 
uniform over A4) in terms of spectral gaps. 

Theorem 4. Let A be an eigenvalue of —Am of multiplicity m, more precisely, 

Afc-i < Aft = A = Afc+ m _i < Xk+m- 
where Xj — Xj(M). Let 5\ = min{l,Afc — Afc_i,Afc_|_ m — Xk+ m -i} and assume 
that p + e/p < CjJ[ k 8\. Let Uk, ■ ■ ■ ,Uk+ m -i be orthonormal eigenvectors of — Ar 
corresponding to eigenvalues Afc(r), . . . , Afc +m _i(T). 

Then there exist orthonormal eigenf unctions gk, ■ ■ ■ , 9k+m—l °/ — A a/ correspond- 
ing to the eigenvalue X and such that 

(7.6) \\uj-P gi \\ 3 <CM,kS?(p + e/p) 

and 

(7-7) |b - L Uj \\ 2 L2 < C M , k 5f(p + e/p) 

for all j — k, . . . ,k + m — 1 . 

Proof. By Theorem [TJ the constant Cm.u in the bound for p + e/p can be chosen 
so that |Aj(r) — Aj(Ajf)| < jS\ for all j < k + m. For every j = k, . . . , k + m — 1, 
apply the second part of Lemma 17.41 with j in place of k, A' = Aj(r) in place of A, 
a = 2\X' - A| and j3 = 7 = |5 A . We have 

X-Sx<X'-j<X<X' + a<X' + (3<X + S\, 
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therefore the assumptions of Lemma 17.41 are satisfied and 

H(x'-y,x'+ a ](M) = H {X }{M) = span{/ fc , . . . , f k + m -i}- 
Thus Lemma IT41 yields that 

(7.8) \\Iu j -g j \\ 2 L2 <C\\'-\\5; 1 + C M>k 5- 2 (p + e/p)<C M ,kS; 2 (p + e/p) 

where (jj = P{\ylv,j. Here the second inequality follows from the fact that |A'— A| < 
CM,k{p + e/p) by Theorem Q] 

By Lcmma l6.2f 1). / is almost isometric (up to an error term Cm,>cP) on the linear 
span of Uk, ■ • ■ , Uk+m-i- This and (|7.8[) imply that the functions . . . , cjk+m-i are 
almost orthonormal up to CM,kS^ 2 (p + e/p). Let {ffj}^™ -1 be the Gram-Schmidt 
orthogonalization of {ffj}^™ 1 , then 

\\9j-9j\\v> <C M ,kS- 2 {p + e/p). 
This and (f778|) imply (|7J| . Now ([7T6]) follows from (|7J| and Lemma [Ol □ 

Note that the functions <?&,... , fffc+m-i constructed in Theorem |4] depend on L 
in rather unpredictable way. The theorem only implies that the subspace gener- 
ated by Iuk, . . . , Iuk+ m -i converges to H^ X }(M) as p + e/p — > 0. A fixed basis 
/ fc , . . . , fk+m-i of H{ A }(M) is approximated by vectors Iu k , ... , Iuk+ m -i trans- 
formed by an m x m orthogonal matrix. 

In the case of multiplicity m = 1, the eigenfunction g k is unique (up to a sign) 
and therefore Theorem |4] implies Theorem [2l 

8. Volume approximation 

This section supplements the main results of the paper. Here we consider various 
aspects of volume approximation in the sense Definition 11.11 

Lemma 8.1 (Marriage lemma for measures) . Let X C M be a finite set. A measure 
p on X is an e- approximation for vol (in the sense of DeHnition ll.l]) if and only 
if vol(M) = p,(X) and for every Y C X one has p(Y) < vo\(U e (Y)). By U £ we 
denote the e -neighborhood of a set. 

Proof. The proof is similar to that of Hall's Lemma for bipartite graphs. The "only 
if" implication trivially follows from the definition. 

We prove the "if" part by induction in N = \X\. To carry on the induction, 
we are proving a more general assertion where M is not necessarily a manifold 
but a metric measure space where the volume of a ball is positive and depends 
continuously on its radius. (In particular, this implies that every sphere has zero 
measure.) Note that the requirement of Definition 11.11 that X is an e-net follows 
from the assumption p(Y) < vo\(U e (Y)) applied to Y = X and the fact that 
vol(Af) = p(X). 

Let X — {xi}fL 1 . The case N = 1 is trivial. Suppose that N > 1 and the 
assertion holds for every metric measure space M' (with the above property) and 
every subset X' C M' of cardinality less than N. We construct a family {Vi(t)}^ 1; 
t £ [0, T], of coverings of M by measurable sets Vi(t) such that 

(1) Vi(0) = B e ( Xi ), and Vi{t) C B e ( Xi ) for all t; 

(2) the sets Vj(T) are disjoint; 

(3) For any set I C {1, ... , N}, the volume of the set \J ieI Vi{t) depends con- 
tinuously on t. 
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Informally, to construct this family, we continuously remove from each Vi some 
pieces of V 3 ■ , j > i. Formally, we set T — N and sequentially define the family for 
t G [0, 1], t G [1,2], . . . , t G [N — 1, N], in such a way that only Vi(t) changes on the 
interval [i — Assuming that the family is already defined for t = i — 1 , we set 

V i {i-t') = V?UB eV {x i ) 

for all t' G [0,1], where V® is the set of points in the set Vi(i— 1) that do not belong 
to any of the sets Vj(i — 1), j ^ i. 

If vo\(Vi(T)) = ^ for all i, then the sets V* = V l {T) satisfy Definition O 
Otherwise consider a maximal interval [0,fo] such that 

(8.i) voi(U iei y i(i) )>5: je/Mi 

for every set / C {1, ...,N}. By continuity, such a to exists and the inequality 
(|8.ip turns into equality for t — to and some set I = Iq C {1, . . . , TV}. (Note that 
(|8.1j) is always satisfied for I = {1, ... , N} since the sets Vi(t) cover M for every i.) 

Let M' = \J ielQ Vi(to) and M" = M \ M' . We apply the induction hypothesis 
to the spaces M' and M" with respective sets X' = {xi} i( zi a and X" = {xi} i( £i , 
equipped with the restrictions of vol and [i. For M' and X', the assumption that 
vo\{M' n U S (Y)) > n(Y) for all Y <Z X' trivially follows from pTT]> . For M" and 
X", we verify the assumption by contradiction. Suppose that 

vol(M" nU E (Y)) < n(Y) 

for some Y C X" . Let Y" = {xi]i e j where J C {1, . . . , N} \ lo- Consider the set 
/ = J U Io ■ For this set we have 

vol fij Fi(t )) < vol(M') + vol(M" n U £ (Y)) < vol(M') + /i(y). 

By the choice of to and Io , we have 

vol(M') = vol (|J V i (t j)='£ Mi 

and /i(y) = Mi by definition. Thus 

vo1 (U ieJ < vol ( A/ ') + ^ = E J6/ 

contrary to (|8.ip . This contradiction proves that M" and X" satisfy the induction 
hypothesis. Thus X' and X" (equipped with the restrictions of //) e-approximate 
M' and M" (equipped with the restriction of vol). Hence {X,fi) e-approximates 
(M,vol). □ 

Recall that the Prokhorov distance |15) 7r(/i, z^) between two finite Borel measures 
/i and v on M is the infimum of all r > such that 

< v{U T {A)) + r and < (J,(U r (A)) + r 

for every measurable set A C M. It is well-known that weak convergence of mea- 
sures is equivalent to convergence with respect to the Prokhorov distance. Let us 
introduce a similar distance 7To(m, v), which makes sense only if /x(M) = u(M). It 
is defined as follows: 7To(m, v ) is the infimum of all r > such that 

< v(U r (A)) and < ^(C/ r (A)) 

for every measurable set A C M. Clearly 7r(/x, ^) < 7r (^, f). 
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Unlike Prokhorov distance, ttq is hardly useful for general metric measure spaces. 
However, in our situation M is a Riemannian manifold and one the measures is its 
Riemannian volume vol. In this case it is more convenient to work with no and it 
defines the same notion of convergence to vol. Indeed, the following lemma holds. 

Lemma 8.2. Let /i be a Borel measure on M such that fi(M) = vol(JVT). Then 



Proof. As already mentioned, the first inequality trivially follows from the defini- 
tions. To prove the second one, let r > be such that 7r(/i, vol) < r. We are to 
prove that, for a suitable r\ = CmT 1 > r an d every measurable A c M, one has 
vol(U ri (A)) > fi(A) and /j,(U ri (A)) > vol(A). If U ri {A) = M, these inequalities 
follow from the assumption that n(M) = vol(M). Suppose that U ri (A) ^ M and 
choose a point p £ M \ U ri (A). Let q be a point nearest to p in the closure of A. 
Connect p to q by a minimizing geodesic and let p\ be a point on this geodesic such 
that \p\q\ = (r\ +r)/2. The triangle inequality easily implies that the ball B r2 (pi) 
of radius r 2 = (r*i — r)/2 is contained in the set U ri (A) \ U r (A). Therefore 



Since r > 7r(/i, vol), we have vo\(U r (A)) > fi(A) — r. Assuming that r is suffi- 
ciently small, we have vol(£? r2 (pi)) ~ v n f2 > r since r 2 > C^vi^ 1 /™. Therefore 
vo\(U ri {A)) > n(A). To prove that fj,(U ri (A)) > vol(A), apply the same argument 
to the set M \ U ri (A) in place of A. This yields that 



where the second inequality follows from the fact that U ri (M \ U ri (A)) C M \ A. 
Since /i(M) = vol(Af), this implies that fi(U ri (A)) > vol(A). Hence 7Tq(^, vol) < r\ 



Proposition 8.3. Let X C M be a finite set and /i a measure supported on X with 
fJ>(M) = vol(Af). Then (X,fx) e- approximates (JW, vol) if and only ifir^lp, vol) < e. 

Remark 8.4. This proposition together with Lemma 18.21 implies that if (X,[i) e- 
approximates (M, vol) then 7r(/x, vol) < e and, conversely, if 7r(//, vol) < ce" (where 
c = c(M)) then (X,fi) e-approximates (M, vol). 

Proof of Provosition \8.3[ First assume that 7To(/x, vol) < e. Then, by the definition 
of 7ro, we have vol(U e (Y)) > [i(Y) for every Y c X. (Here we use the fact that 
the boundary of U e (Y) has zero volume.) This and Lemma \8 . 1 1 imply that (X,/j,) 
e-approximates (M, vol). 

Now assume that (X, fj,) e-approximates (M, vol) and let A C M be a measurable 
set. It suffices to prove that 



7r(jU,vol) < 7To(/x, vol) < Cx7r(/x, vol) 1 



vol(U ri (A)) > vo\{U r (A)) + vo\(B r2 ( Pl )). 



H(M \ U ri (A)) < vol(U ri (M \ U ri {A))) < vol(M \ A), 



and the lemma follows. 



□ 



(8.2) 



vol(U e (A)) > fi{A) 



and 




fi(U e (A)) > vol(A). 



vo\{U e (A)) > vo\{U e {A n X)) > (j,(A nX) = fj,(A) 
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where the second inequality follows from Lemma 18.11 To prove (|8.3|) , consider the 
set Y = X n U £ {A). By Lemma IO we have 

(*(X \Y)< vol(U E (X \ Y)) < vol(M \ A) 

where the second inequality follows from the fact that U S (X \ Y) C M \ A. Since 
fi(X) = vol(Af), this implies that fi(Y) > vol(A) and ((83)) follows. □ 

Computing the weights. We conclude this section by discussing how weights /Xj 
can be computed, for a given an e-net X — {x{\ C M . There are several natural 
ways to associate a partition {V{\ as in Definition 11.11 to the e-net. One is to let 
{Vi} be the Voronoi decomposition of M determined by X. Another possibility is 
to define Vi = B e (xi) \ U,<, ^ e {xj). However actual computation of the weights 
= vol(Vi) is these constructions may be complicated. 

A more practical approach could be the following. First split M into small 
subsets (of diameter at most e' < e) whose volumes are easy to compute. To each 
of these subsets, associate one of the nearby points from X. Then the weight \Xi can 
be defined as the sum of volumes of the subsets associated to the point Xi. These 
weights define an (e + e')-approximation of volume. 

For example, let M c W 1 be a bounded region (rather than a closed manifold). 
Then M (except a small neighborhood of the boundary) can be divided into small 
coordinate cubes. To each cube one could associate the point of X nearest to the 
cube's center. The resulting weights are roughly equal to the volumes of Voronoi 
regions but are easier to compute. 

Remark 8.5. It is an interesting problem how to derive the weights from the dis- 
tance matrix of X without referring to the manifold M. Ideally, one wants a nice, 
symmetric formula for fa in terms of the distance matrix. We were not able to 
come up with such a formula. However there is a straightforward algorithm based 
on the property that a Riemannian metric is locally almost Euclidean. 

Let r £ (Ce,p), Kr 2 <C 1. Then r-balls in M are bi-Lipschitz close to the 
r-ball in W 1 . Moreover for each point Xi £ M one can construct a bi-Lipschitz 
almost isometry <p: B r (xi) — ¥ K n from distance functions of points xj £ B r (xi). 
For example, a function of the form x t-> d(x, Xj) 2 — d(x, Xi) 2 is close to a linear one 
in geodesic normal coordinates centered at Xi . Using such functions as coordinates 
and post-composing with a suitable linear transformation of W 1 one gets a desired 
almost isometric map (p. The image cp(X D B r (X)) C K n is easy to compute, and 
then the problem essentially reduces to computation of volumes of Voronoi regions 
(or differences of balls) in R™. 
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